The electron-phonon coupling is large for localized states 
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From density functional calculations, we show that localized states stemming from defects or 
topological disorder exhibit an anomalously large electron-phonon coupling. We provide a simple 
analysis to explain the observation and perform a detailed study on an interesting system: amor- 
phous silicon. We compute first principles deformation potentials (by computing the sensitivity 
of specific electronic eigenstates to individual classical normal modes of vibration). We also probe 
thermal fiuctuations in electronic eigenvalues by first principles thermal simulation. We find a strong 
correlation between a static property of the network [localization, as gauged by inverse participation 
ratio (IPR)] and a dynamical property (the amplitude of thermal fluctuations of electron energy 
eigenvalues) for localized electron states. In particular, both the electron-phonon coupling and the 
variance of energy eigenvalues are proportional to the IPR of the localized state. We compare the 
results for amorphous Si to photoemission experiments. While the computations are carried out for 
silicon, very similar effects have been seen in other systems with disorder. 
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I. INTRODUCTION 

Electron states of finite spatial extent, so called "local- 
ized" states are ubiquitous in nature. Localized molecu- 
lar orbitals occur in molecules and in solid state systems 
with disorder (this could take the form of defects, such 
as a dangling bond state or more subtly from topologi- 
cal and or chemical disorder in amorphous materials or 
glasses). Such localized states are also known in poly- 
mers. Surface states, as another form of defect, also are 
often spatially compact. The study of localization in its 
own right has been an important and active subfield of 
condensed matter theory since the fifties. 

Another physical quantity of key importance is the 
electron-phonon (e-p) coupling, the interaction term con- 
necting the electronic and lattice systems. Perhaps most 
spectacularly, the e-p coupling is the origin of super- 
conductivity as expressed in BCS theory^. Phillips^i has 
shown that large e-p couplings in the cuprate super- 
conductors can lead to a successful model of high Tc 
superconductivity"^ within the framework of conventional 
BCS superconductivity. The e-p coupling is also the me- 
diator of all light-induced structural changes in materials. 
In amorphous silicon the greatest outstanding problem of 
the material, the Staebler-Wronski effect"', depends criti- 
cally upon the electron-lattice interaction. A zoo of anal- 
ogous effects is studied in glasses; perhaps the most im- 
portant example is reversible photo-amorphization and 
photo-crystallization used in the GeSbTe phase-change 
materials used in current writable CD and DVD technol- 
ogy. 

Previous thermal simulations with Bohn-Oppenheimer 
dynamics have indicated that there exists a large 
electron-phonon coupling for the localized states in the 
band tails and in the optical gap^-^. Earlier works on 
chalcogenide glasses by Cobb and Drabold^ have em- 
phasized a strong correlation between the thermal fluc- 
tuations as gauged by root mean square (RMS) vari- 



ation in the LDA eigenvalues and wave function lo- 
calization of a gap or tail state (measured by inverse 
participation ratio^, a simple measure of localization). 
Drabold and Fedders^ have also shown that localized 
eigenvectors may fluctuate dramatically even at room 
temperature. Recently, Li and Drabold relaxed the adi- 
abatic (Born-Oppenheimer) approximation to track the 
time-development of electron packets scattered by lat- 
tice vibrations-. In this paper, we examine the electron 
phonon coupling and provide a heuristic analysis of the 
e-p coupling for localized electron states. We explore the 
e-p coupling in some detail for a particular model sys- 
tem (amorphous silicon) which provides us with a con- 
venient variety of localized, partly localized "bandtail" 
and extended states. We compute deformation poten- 
tial (which measures the response of a selected electron 
state to a particular phonon), and also track thermally- 
induced fluctuations of electronic eigenvalues. We find 
that localized states always exhibit a large e-p coupling. 
Our computations are carried out using a first princi- 
ples molecular dynamics code Siesta, and the eigenval- 
ues and states that we study are from the Kohn-Sham 
equations with a rich local orbital basis. A rationale for 
the study of the Kohn-Sham states is given elsewhereiSi. 
We emphasize that the results that we give are qualita- 
tively general - not just an artifact of studying a disor- 
dered phase of silicon (we have, for example, seen exactly 
the same effects in various binary glasses which exhibit 
very different topological and chemical disorder). 



II. THEORY 

To establish a connection between electron-phonon 
coupling and wave function localization for the electrons, 
we consider an electronic eigenvalue A„ near the band 
gap of a-Si. The sensitivity of A„ due to an arbitrary 
small displacement of an atom (possibly thermally in- 
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duced) can be estimated using the Hellmann-Feynman 
theorenJi, 

Here we have assumed that the basis functions are fixed 
and \ipn) are the eigenvectors of the Hamihonian H. For 
small lattice distortion {(5Rq}, the corresponding change 
in 6Xn is, 
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where N is the total number of atoms in the model. If 
the displacement (5Rc((t) arises from classical vibrations, 
one can writei^, 



3N 



5R„(t) = ^A(T , Uj) COs(w t + (f)^) Xa{^), (2) 

where u indexes the normal mode frequencies, A(T, lu) is 
the temperature dependent amplitude of the mode with 
frequency oj, (j)^ is an arbitrary phase, Xq(^) is a normal 



mode with frequency u and vibrational displacement in- 
dex a. Using the temperature dependent squared ampli- 
tude — SkeT/Maj^, the trajectory (long time) average 
of SXf-^ can be written (using Eq.nand[21) as. 



where the electron-phonon coupling ^n(^) is given by, 
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One can infer from Eq.El that thermally induced fluc- 
tuation in the energy eigenvalues is a consequence of 
electron-phonon coupling. Note that for a given elec- 
tronic eigenvalue, the contribution to the coupling comes 
from the entire vibrational spectrum involving all the 
atoms in the systems. Since the normalized eigenstate 
can be written as \ipn) = X)i N) where \i) are the ba- 
sis orbitals, it follows from Eq.^that, 
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The first term in the second line of Eq.|31 is positive 
definite (diagonal) while the second one, the off-diagonal 
term (indicated by the prime), is not of a single sign. 
In the event that only a few a„i dominate (the case for 
localized states), then the leading contribution to the 
electron-phonon coupling originates largely from the di- 
agonal term. The addition of a large number of terms of 
mixed sign and small magnitude leads to cancellations in 
the off-diagonal term leaving behind a small contribution 
to electron-phonon coupling. By comparing to direct cal- 
culations with the full Eq.|31 we show that dropping the 
second term appears to be reasonable for well-localized 
electron states. The approximate "diagonal" electron- 
phonon coupling can be written as. 
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where Nf, is the number of basis orbitals and qni — 
fln.iP is the charge sitting on the ith orbital for a given 
normalized eigenstate {ipn}- The degree of wave function 
localization can be measured by defining inverse partici- 
pation ratio X for the eigenstates |V'n), 
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X(n) = qni 
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Equation^leads to an approximate but analytic connec- 
tion between X and electron-phonon coupling. Since X is 
large for locahzed states, one expects S„(w) (and there- 
fore (SXn)) to be large for a locahzed state. If we fur- 



ther assume that 7^(ijj,i) = X]Q=i[(*I^H/9RQ|i)]^Xa('^) 
is weakly dependent upon site/orbital index i, then 



(8) 



where /(w) is defined from 7^ and Eq.|Sl In this "sepa- 
rable" approximation, it is also the case that (SX^) oc Z„. 
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FIG. 1: Electron-phonon coupling surface plot for a 216-atom 
model of a-Si. The absolute value of electron-phonon coupling 
|H| (cf. Eq.0J is plotted as a function of phonon frequency uo 
and energy eigenvalues near the gap. The largest value of |H| 
in the plot corresponds to the eigenvalue for HOMO, which 
is the most localized state in the spectrum. 



FIG. 2: (Color online) Mean square fluctuations of electronic 
eigenvalues versus inverse participation ratio plot at differ- 
ent temperature. The fluctuations at temperature 150K and 
300K are found to be linearly correlated with the participation 
ratio for the corresponding eigenstates as predicted in section 
II. The correlation coefficient (r) for different temperature is 
indicated in the plot. 



III. METHOD 

The model of a-Si we have used in our calculations 
was generated by Barkema and Mousseaui^ using an 
improved version of the Wooten, Winer and Weaire 
(WWW) algorithmic. The details of the construction 
was reported in Ref.^3. The model consists of 216 atoms 
of Si packed inside a cubic box of length 16.282A and has 
two 3-fold coordinated atoms. The average bond angle is 
109.5° with a root mean square deviation of 11.0°. The 
density functional calculations were performed within the 
local density approximation (LDA) using the first prin- 
ciples code SlESTAiSiiSiii. We have used a non self- 
consistent version of density functional theory based on 
the linearization of the Kohn-Sham equation by Harris 
functional approximationM, along with the parameteriza- 
tion of Perdew and Zunger— for the exchange-correlation 
functional. The choice of an appropriate basis is found 
be very important and has been discussed at length in a 
recent communicationifi. While the minimal basis con- 
sisting of one s and three p electrons can adequately de- 
scribe the electronic structure of amorphous silicon in 
general, there is some concern about the applicability of 
these minimal basis in describing deeply localized and low 
lying excited states in the conduction bands accurately. 
We have therefore employed a larger single-C basis with 
polarization (d) orbitals (SZP)^°'^^ in the present work. 



Throughout the calculation we have used only the F point 
to sample the Brillouin zone. 

Starting with a fully relaxed configuration, we con- 
struct the dynamical matrix elements by successively dis- 
placing each atom in the supercell along three orthogonal 
directions (x, y and z) by O.OlA and computing the forces 
for each configuration. Within the harmonic approxima- 
tion, the spring constant associated with each atom and 
direction can be written as a second derivative of the to- 
tal energy with respect to the displacement of the atom in 
that direction. We have checked the convergence of both 
the matrix elements by using a different set of values for 
atomic displacement and used a value of O.OlA in our 
calculations. 9A/9Rq. was obtained by finite differences 
from the dynamical matrix calculations. The calculation 
does not need any extra effort beyond that of forming 
the dynamical matrix. 

To explore the validity of our analysis and to eluci- 
date the connection between the localization (IPR) (I) 
of electronic eigenstates and fluctuation of the conjugate 
eigenvalues, we performed thermal MD simulations at 
constant temperatures using a Nose-Hoover thermostat. 
The simulations were performed at temperatures 150K, 
BOOK, 500K and TOOK with a time step of 2.5fs for a to- 
tal period of 2.5ps. For a given temperature, the mean 
square fluctuations were computed by tracking the eigen- 
values at each time step and averaging over the total time 
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FIG. 3: The average electronic density of band tails states for 
four different temperature T = 150K, 300K, 500K and TOOK. 
Note that the conduction band tails (right) near the Fermi 
level (which is between the two tails) are more sensitive to 
thermal disorder than the valence band tails (left) providing 
a qualitative agreement with experimental result in Ref. 25. 

of simulation excepting the first few hundred time steps 
to ensure equihbration. The mean square fluctuation (TV) 
for an energy eigenvalue An(t) is defined as : 

7^„ = ((A„(t) - (A„))2), (9) 

where ( ) denotes the average over time. We study the 
fluctuations of {An(t)} by plotting against time at a given 
temperature and compare it with the T obtained for the 
corresponding eigenvalues. For illustrations of such adi- 
abatic evolution of Kohn-Sham eigenvalues, see Ref. 0. 

IV. RESULTS 

In Fig.^ we have plotted the electron-phonon coupling 
for the states near the band gap obtained directly from 
Eq.0] It is clear from the figure that the e-p coupling is 
large only in the vicinity of conduction and valence band 
tails. The largest e-p coupling in the plot corresponds to 
the highest occupied molecular orbital (HOMO) in the 
optical frequency regime around 415 cm"'^. The low- 
est occupied molecular orbital (LUMO) also has a large 
feature around the same frequency. A MuUiken charge 
analysis and inverse participation ratio calculation of the 
electronic eigenfunctions have shown that both these two 
states - the HOMO and LUMO are highly localized and 



are centered around the dangling bonds present in the 
model. On moving further from the band tails in either 
direction along the energy axis, the e-p coupling drops 
quickly and the surface becomes featureless for a given 
eigenvalue. This behavior of e-p coupling can be under- 
stood from the arguments presented in section II where 
we have shown that the e-p coupling for localized states 
is directly proportional to the inverse participation ra- 
tio. For a localized state, therefore, the large value of 
electron-phonon coupling can be attributed to the large 
value of inverse participation ratio associated with that 
state. Since HOMO and LUMO are the two most lo- 
calized states in the spectrum, the e-p coupling is large 
for these states and as we move toward the tail states, 
the coupling decreases. It is important to note that the 
plot in the Fig.^ has been obtained from Eq.^ with- 
out making any approximation and is exact inasmuch as 
the matrix elements obtained from the density functional 
Hamiltonian are correct. This observation supports our 
assumptions that the dominant contribution to e-p cou- 
pling comes from the diagonal term in Eq.|31 and that 
7^(a;,i) is weakly dependent upon site/orbital index i 
and also indicates from direct simulation there exists a 
linear relationship between mean square fluctuation of 
electronic eigenvalues and the corresponding inverse par- 
ticipation ratio for localized states. 

In order to justify our arguments further presented in 
section II, we now give a look at the mean square fluctua- 
tion of energy eigenvalues. As outlined in section HI, we 
have computed the mean square fluctuations at four dif- 
ferent temperature (150K, 300K, 500K and TOOK) from 
MD runs over a period of 2.5ps and plotted in the Fig.|5| 
The fluctuation obtained this way provides a dynamical 
characteristic of the band tails states and is compared 
with a static property, the inverse participation ratio of 
the same states. A simple linear flt reveals a strong cor- 
relation between the eigenvalue fluctuation and the cor- 
responding inverse participation ratio for the states. The 
correlation is found to be as high as « 0.95 for T=150K 
and 300K and falls to w 0.8 at high temperature. The 
value of the correlation coefficient for different tempera- 
ture is indicated in the Fig.[21 Once again, we see that 
the result is in accordance with our prediction in section 
II and provides a simple physical picture for having a 
large electron-phonon coupling for the localized states. 

In Fig.O we have plotted the time averaged electronic 
density of states for four different temperature in order 
to study the effect of thermal disorder on the tail states. 
It is quite clear from the figure that the effect of ther- 
mal broadening is quite significant on both sides of the 
gap. Photoelectron spectroscopic studies on a-Si:H by 
Aljishi et al.— have shown that the conduction tail is 
indeed more susceptible to thermal disorder than the 
valence tail. The temperature dependence can be con- 
veniently expressed by introducing a characteristic en- 
ergy Eq and fitting the electronic density of states to 
p(E) w exp(|E - Ef |/Eo(T)). Aljishi et al. expressed the 
temperature dependence of the tail states by the slope 
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of the Eo(T) vs. T plot and obtained a smaller value 
for the conduction band tail. We have observed a qual- 
itative agreement of our results with experiment. The 
key observation that one should note from Fig.O is the 
following: the shape of the tail in the conduction band 
rapidly changes as the temperature rises from 150K to 
700K. The corresponding change in the valence tail for 
the same range of energy (0.4eV) is however much less 
and is rather smooth compared to the conduction tail. 
Since the localized defect states (coming from the two 
dangling bonds) have been removed before plotting, this 
observation qualitatively suggests that the conduction 
tail states are more susceptible to lattice motion. It is 
tempting to attempt to estimate decay parameters for a 
direct comparison to experiment—, but the sparse sam- 
pling of tail states for this 216-atom model makes this a 
dangerous exercise. The basic features do appear to be 
represented in our study, however. 

V. CONCLUSION 

Using accurate methods and a reasonable model of 
a-Si, we showed that there is 1) a large e-p coupling 



for localized states, 2) a significant correlation between 
thermal fluctuation of electron energy eigenvalues con- 
jugate to localized states and the IPR of the model at 
T=0, 3) We find qualitative agreement with photoemis- 
sion experiments^^ , 4) we provide a simple analytic argu- 
ment for the origin of these effects. Identical experience 
with models of other amorphous materials has convinced 
us that the results are correct in at least a qualitative way 
for binary glasses and amorphous materials, and perhaps 
other systems beside. 
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